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ABSTRACT 

The present paper considers optimal means of characterizing the distri- 
bution of product energy states resulting from reactive collisions of molecules 
with restricted distributions of initial states, and vice versa, iV£., character- 
izing the particular reactant state distribution which yields a given set of 
product states, at a specified total energy E . The S-matrix, or reaction 
probability matrix P(E) , "global" in nature, contains much more detail than 
necessary to reproduce the results of any single specific experiment or computer- 
simulation thereof (via classical mechanical trajectory calculations). More- 
over, since reactant and/or product state resolution is always experimentally 
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limited (to a greater or lesser degree), data are necessarily coarse- 
grained accordingly. Many quantal features are thereby lost and the 
results are often at a level appropriate for comparison with classical 
calculations (e.g .,in the form of low-resolution contour maps of energy 
disposal) . Such contour plots of the yield function Y or the averaged 
transition probability to (the "poor-man’s" ^“matrix) nevertheless con- 
tain the essence of the dynamical results. It is suggested to represent 

l 

the energy-dependence of global-type results in the form of square-faced 
bar plots, and of data for specific-type experiments (or computer simula- 
tions) as triangular-faced prismatic plots (contour maps vs. E ). . The 
essential parameters defining the internal state distribution are isolated 
and the information content 1(E) of such a distribution (for a micro- 
canonical ensemble) is put on a quantitative basis. The relationship 
between the information content, the surprisal, and the entropy of the 
continuous distribution is established. The concept of an "entropy 
deficiency" AS* , which characterizes the specificity of product state 
formation, is suggested as a useful measure of the deviance from 
statistical ("phase— space dominated") behavior. The degradation of 
Information by experimental averaging is considered, leading to bounds 
on the entropy deficieheyc 
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1. Introduction 

Important information on the dynamics of reactive molecular collisions 
is derived from a knowledge of the distribution of internal energy states of 

the products or the relative reactivity of reactants in different excited 

1 21 

internal states . 

In the present paper, for simplicity we restrict our attention to the 
simplest case of a three-center, atom-transfer reaction of the type 
A + BC AB* + C, where the dagger denotes internal excitation of the product 
diatomic. An example is the hydrogen- abstraction reaction: 

Cl + HI (v,J) + HC)l(v\J') + 1, (1) 

one of a number of three-center reactions which have been well-studied via 

the infrared chemiluminescence technique by Polanyi and his followers 1-5 . 

Ordinarily the initial internal state (v,J) distribution has been 

essentially Boltzmann, fixed by the temperature of the reactant molecules BC. 

The experimental measurements yield the internal state distribution (v',J') 

f 

of the product molecules AB , often found to be non-Boltzmann in character. 

In many cases extensive population inversion occurs, thereby providing the 
basis of chemical laser action . 

With the exception of the hydrides, the resolution of individual 
internal quantum states of the products has not yet been feasible. For most 
practical purposes, when the density of internal states is high, it is 
therefore convenient to regard the internal energies as continuous variables, 
say Ey, and E R , , representing product vibrational and rotational energies, 
respectively. Even when individual states have been resolved, it is often 
useful to "smear out" the discrete distribution to convey a qualitative 
picture of the overall experimental results. Moreover, with the wide 
utilization of classical mechanical trajectory calculations and Monte Carlo 
evaluation of reactive scattering cross sections it behooves us to adapt 
to classical, continuous probability density functions, histograms and "bins", 
rather than to attempt to "quantize" artificially the classically calculated 
results. In sections 2-4 and 9 we will concern ourselves mainly with coarse- 
grained product (and reactant) energy distributions. 

Given the internal energy distribution of product (AB^) molecules 
one can immediately calculate the distribution of relative translational 



energy „ E^, , 
relation : 


o£ the recoiling products , making us® of the energy conservation 





= E T * E V + E R " 
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energy of the A ■> BC system, E^, the vibrational and rotational energies , Ey 

and E^ s of the reactant molecule BC, and the zero-point to zero-point 

reaction exoergicity, “AE q . The translational distribution of products can 

be measured directly , in the course of molecular beam scattering experiments , 

7=9 

via the velocity analysis technique . The crossed molecular beam method 

10 

has also yielded more direct information on product state distributions 

For a number of reactions it has been possible to carry out rotational state 
11 12 

analysis, s and, more recently, vibrational state analysis has been 

accomplished 1 s . The chemical laser technique, as pioneered by Pimentel 

and co-workers, has provided independent data on relative populations of 

1 15 

vibrational states of product molecules, s which, ar® directly comparable 
with the ir chemi luminescence results . Xt is clear that the question of the 
energy disposal in elementary exosrgic reactions is ©:? great imports, o© in 
the field of modem chemical kinetics. 

Alternative information on microscopic chemical reaction dynamics 
is available from experiments involving reactants in excited internal states. 
Thus far such studies have been concerned with endoergic reactions or those 
with an energy barrier to be surmounted. The influence of reactant internal 
energy ^ s + E^, controlled by the temperature of the BC molecule, at 

constant E v , has been studied in one such case by Schmeltekopf et al . and 
the separate effect upon the reaction cross section of E.. at fixed E_ and E r 

18 V K 19 k 

has been investigated by Chupka e t al . , and by Brooks and co-workers 
The possibility of utilizing judicious vibrational excitation of reactants 

for selective influence of the product molecular configuration has already 

20 

been considered , and selective excitation by lasers is now a rapidly 

21 22 
developing field with many practical implications. 

The considerable experimental activity in this field has stimulated 

23 24 

much theoretical interest 4 , Anlauf et al . pointed out the utility of 
microscopic reversibility to gain information m th© influence of reactant 



internal state distributions for the exoergic reaction. Kinsey showed how 

such considerations could be most fully exploited when only limited data of 

26 

various types on 1 the Exoergic reaction are available. Marcus further 

investigated these same' questions and showed the relationship to the quasi- 

23 27 

equilibrium approximation. ■ The simplistic optical model analysis 5 has 

also been used in this “area, interrelating the forward and reverse reaction 

28- 30 

cross sections . 'The statistical theory has been found to be a useful 

diagnostic tool,, since deviations from its predicted {equilibrium micro- 
canonical) distribution of product states are a measure of the specific, non- 
equilibrium nature of the reactive collisions. 

The ultimate theoretical description of the dynamics is, of course, 
a proper dynamical theory. At the present stage there are still formidable 
obstacles in the path of a full quantal solution of the reactive scattering 
problem. 31 1 Thus we turn to classical mechanical, numerical trajectory 
calculations i ’ ’ which attempt to stimulate experiment assuming a 

"realistic", adiabatic potential surface and neglect all quantal interference 
effects. In what follows we shall assume that at least fragmentary, low- 
resolution data or computer- simulated experimental results are available, and 
our goal is to extract as much information as possible from the available 
"data", whether experimental or calculated. 

One of the aims of the present paper is to develop optimal means 
of characterizing the distribution of product states (in a given experiment, 
at a given total energy) , as governed by the reactant state distribution (and 
vice versa, i.e . , the reactant state distribution yielding a given set of 
product states), in terms of any two of the three energy variables E^,, , Ey, , 
E_, . Of course, only two of these are independent quantities (cf.Eq. 2). 

It will often be convenient to work with fractional quantities f x , defined 

as the fraction of the total available energy in the mode X (where X may 
be T,V,R or the same symbols primed). Throughout the paper attention is 

restricted to the energy range below the threshold for collision-induced 

dissociation. 

Before proceeding further, it is necessary to distinguish between 
two classes of data, i.e., "global" and "specific". For a complete statement 
of the dynamics of the reaction, one requires the transition probabilities 
from all possible initial states to all possible final states. Such a global 
representation requires a transition probability matrix, ]?, whose elements 
are mod-squared elements of the 5-matrix, such as would be provided by a full 



quantal solution to the scattering problem. Neither the experimenter nor 
the classical mechanical computer are as ambitious as the quantum mechanic, 
however. They are content to describe an incomplete but specific experiment, 
i . e . , a distribution of product states (or product energies) for a given 
set of initial conditions. 

One can further distinguish (within the class of "specific" 
representations) between a detailed study in which the product distribution 
is determined for given, specific initial reactant states, and an " inclusive" 
study which yields the product state distribution averaged over the reactant 
state distribution. Assuming perfect resolution of internal states, a 
global description provides the whole P- matrix. A detailed, specific study 
(later referred to as "elusive") characterizes a column (or a row) of the 
Py matrix, while an inclusive, specific study provides only an "average" 
over a row or column. 

Sec. 2 considers the presentation of the results of a study of 
the global dynamics in the form of square-faced "bar" plots; such a plot 
provides low- resolution probability contour maps (the poor-man's P^-matrix) . 
Sec 3 discusses the representation of "specific" results in the form of 
a triangular- faced "prismatic" plot. The essential parameters defining 
a given internal state distribution are isolated and discussed in Sec. 4„ 

The "information content" of a product state distribution is considered in 
Sec. 5. Transition probabilities and the problem of degenerate (or 
experimentally indistinguishable) states are discussed in Sec. 6 (which is 
confined to discrete distributions). Sec. 7 is concerned with the 
information content of the two types of "specific" experiments, as well 
as the "global" type of results. Sec. 8 summarizes the hierarchy of 
possible representations that are available for the distribution of products 
(and influence of the distribution of the reactants). Sec. 9 considers 
the problem of defining the information content of a continuous distribution. 
Sec. 10 makes the connection with entropy and suggests the concept of an 
entropy deficiency characterizing the specificity of the product state 
distribution. Sec. 11 summarizes the important definitions and results. 
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2. SQUARE-FACED BAR PLOTS FOR SUMMARIZING RESULTS OF "GLOBAL" EXPERIMENTS 

One may "collapse" the transition probability matrix P, at a given E, 

to a function of two independent variables E^, E T , (or E i n t < are 

cony lement ary variables). Primes designate products (post-collision variables). 

The fractions fj, £j t are convenient reduced variables; thus f T ■ E-jV^T^max 

and f T , = E T ./( E T .) max » where ( E T ) max » ^T'^max are the maximum conservation- 

allowed values of E t> E^,, respectively (cf. Sec. 1). Thus 1-f^, and 1-f^., are 

the reduced internal energy of reactant and product molecules, respectively. 

25 

Let to(E 1 ,,E T ,) be the average transition probability corresponding to 
reactants with E T in the range between E T and E^ + dE^ forming products with 
E^ # in the range E T , , E,j,, + dE T , (at a given total E) . Then 

»Cb t ,e t ,) - j'j} (r Kn Vi' : % O) • (3) 

Here n,n' denote the internal states of the reactants and products respectively, 

2 

P , a Is , I is the n -*-n' state-to-state reaction probability and the sums 
n' ,n 1 n' ,n 1 r 

are confined to states n,n' in the specified ranges (corresponding to the 
E T ,E T , increments). The concept of averaged transition probability (and its 
symmetry with respect to the interchange of E^ and E T ,) is further discussed 
in Secs. 6 and 9. It is of course clear that the concept of a smooth functional 
dependence of w on the two variables Ej, E^,, implies not only classical mechanics 
but also a sufficiently high density of internal states so that the dependence 
on the variables is smooth and the limit implied by (3) (both numerator and 
denominator are differentials) is well defined. Hie probability ia can be 
presented as a contour map in terms of E^,, E T , (as in Fig. la) or in terms 
of f T , f T , on a unit square (Fig. lb). Of course we must take cognizance 
of the relation between the two 

w (E x )dE x - . u(f^jdf x , (4) 

Note that, in general, io(E t ,E t ,) is the same quantity for both the forward 
and the reverse reaction, and a single contour map is sufficient to characterize 
the global dynamics at a specific total energy E. Stacking together such squares 
at successively greater values of E yields a 'bar plot", whose cross section is 
the f.p, f^, plane and long axis and the E scale. The loci of the maximum of 
(coordinates f^ m ^ , fj^) are plotted vs. E in Fig. lc. 

There are no "global" results available yet, either experimental or 
computer- generated, even at a single E, although a considerable amount of 
fragmentary information from different sources is known. For example, data 
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exist which correspond to vertical or horizontal ’’slabs" of the bar. 

In certain common types of "inclusive" experiments the total 
reaction cross section into all product internal states is measured as 
a function of E^, with reactant internal state distributions essentially 
Boltzmann,, corresponding usually to small £y, f^. Thus tj, is maintained 
at a near-constant value slightly less than unity while E (= E T - AE ) 
varies directly with E^. This measurement of o(E^,) is the most common 
one in conventional beam experiments , for example in studies involving 
ionic reagents, * collisional ionization of neutrals, and, 

recently, for neutral reactants yielding neutral products^*. Such 
experiments amount to summing (or averaging) over all £„, in a "vertical 
slab" of a bar-plot such as Fig. lb. 

There are more refined, and thus only partially-inclusive, 

experiments in which the product translational distribution is measured 

■ 7b 

as a function of E™ (with the reactant internal state distribution held 
constant). Here the "vertical slab" has been analyzed in terms of its 
distribution. 

In another type of experiment the reactant internal energy is varied 

at nearly constant E (or over a small range of E) and the total reaction 

1 ' 7 18 

cross section observed," * s ' but with no translational (or internal) 
energy analysis of the products. These results correspond to summing over all 
fp, for a series of successive vertical slabs at given f^, values. An example 
of such primary data is shown in Fig. 2, constructed from the results of 
Refs. 17 and 43. 


The global results require observations (or computations, at least 
via classical trajectories) at constant E over an entire square grid of 


£p fp, to establish a contour- map at each of Several values of E. Such a 
global bar-plot would have practical indications with regard to 
optimization of reaction yield. The relative importance of E^ vs_. at 

v Cl 

given E can be ascertained by inspection of such^ contour map. E. g. , from 

Fig. la (at E = 10 units) , optimum yield of products for the endoergic 

reaction would be obtained if the reactants had a value of E^^ of 6 units 

(i.e., = 0.4^ thus E^™^ - 4 and = 6 units) . Note that the smallest 

"allowed" value bf’E'is -AE = S units. Information on both forward and 

o 

reverse reactions is contained in the same contour map. The limitation of a 
Square =£ ace d 'bar-plot such as that of Fig. lb is, of course, that distinction 
is made only between translational and internal energies, whereas we know 


that the separate internal modes (rotation and vibration) contribute 
very differently to the reaction probability. Thus we should consider 
the "decomposition" of E^^. into its two "components", Ey and E^, and 
their separate effects, as discussed in the next section. 
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3. TRIANGULAR- FACED PRISMATIC PLOTS FOR SUMMARIZING RESULTS 

OF "SPECIFIC" EXPERIMENTS 

The information content of one type of "specific" experiment, e. g. „ 
on© in which reactant states are specified, is the coarse-grained 
distribution of product states represented as functions of any two 
independent variables of the set {E^, , E y , , E R , } . We may also consider 
the inverse of this, i . e. , the distribution of reactant states, in terms 
of {E t , E yf E r } , which yield a specified distribution of product states 
corresponding to the same total energy E. If this specified product state 
distribution is chosen to coincide with the reactant state distribution of 
the previous experiment, the results of the two experiments are uniquely 
related via the principle of microscopic reversibility. 24 " 26 The properly 
symmetrized measure of reaction probability is either the yield function, 44 

o r 

Y, or the closely- related averaged transition probability, w. 

The results of a "specific" experiment can be conveniently represented 

as a contour map in a plane, with contour lines denoting equiprobable final 

24 

states (or initial states, as the case may be). Of the several possible 
choices of diagrams (see Appendix 1 for details), the equilaterial triangular 
one used by Kinsey seems preferable, since it maintains complete symmetry 
with respect to the three modes T,V,R (or T',V' ,R'). Such a plot is 
illustrated schematically in Fig. 3. The coordinate axes are the fractions f 

x 

of the available energy in mode X, and the vertices designated T, V, R 
represent 100% of the energy in the specified mode, as usual. Of course, any 
two coordinates of the set (f^., fy, f R ) suffice to specify the location of 
a given point. (Note that such a triangular probabi lity contour map is not 
the same as a triangular map of relative rate constants, as used to display 
the chemiluminescence results . ** 24 ’ 48 ) 

The most significant single characteristic feature of such a 
continuous state distribution or contour map of reaction probabilities 
is the location of the maximum probability (and possibly of any subsidiary 
maxima, not present in the example shown) . The coordinates of the maximum 
are designated , fW any two of which fix the "most probable" 

location. The next most significant feature of such a state distribution is 
the "sharpness" of the maximum,. This question is discussed in Sec. 4. 
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Of course the location of the maximum and its intensity (and curvature) 
will depend upon E, but in general one expects that the energy dependence 
of the most-probable fraction f^ m ^ in mode X will be less severe than that 
of , the most-probable energy^'^ itself. 

In order to display the overall energy dependence of the state 
distribution one can array a "stack" of such contour maps along a total 
energy axis, E, perpendicular to the triangular faces, thereby obtaining 
a prism, as shown in Fig. 4a. Each cut through the prismatic diagram 
represents the energy partitioning at the specified, constant total E. 

Fig. 4b summarizes the most significant data, the trajectory of the locus 
of the maximum yield as a function of E, analogous to that of Fig. lc. 

It should be noted that the entire prismatic plot corresponds to 
a certain set of initial conditions specified by fy, f R or, more commonly, 
a particular distribution in fy, f R . We may thus distinguish between two 
types of "specific" experiments. The first is an "inclusive" type, where 
there is a broad distribution of initial reactant conditions ( i. e .. 
distribution in fy, f R ) which give rise to the product distributions in 
fy, , f R , , fp, displayed on each triangle. Conversely, for a given 
specified fy, f R one has data only on the total product yield and so the 
contour map in the triangle summarizes the influence of the initial state 
distribution in f v , f R upon the overall reaction probability. An example 

v K j 

of the former is the usual ir chemiluminescence experiment in which 
Boltzmann distributions of reactants at given temperatures describe the 
initial conditions and then the experimental results consist of detailed 
product state distributions, i .e. . u>(fy, , f RI ) valid for the set {Fy, f R ) 
at some mean value E of the total energy. An example of the inverse is 
an experiment in which the total reaction cross section (irrespective of 
product states) is measured as a function of Ey and/or E R at essentially 
fixed E, as in Refs. 17-18. ' - _ _ 

It is clear that the results of such "inclusive", "specific", 
experiments provide only a limited amount of information on the reverse 
reaction. This can most easily be seen upon consideration of the square- 
faced bar-plots of Sec. 2. An inclusive experiment supplies only the (possibly 
weighted) average (or sum) of entries along a row or column and hence cannot 
characterize the entire square face. 
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As an example, consider a simple case of an exoergic reaction 

(cf. Eq. 1) carried out in a crossed beam experiment, with reactant BC 

in a Boltzmann internal state distribution (characterized by some temperature 

Tint^ 311(1 with 3 g iven relative translational energy distribution, say 

Maxwellian (corresponding to some temperature T ) . Suppose that this, 

inclusive, experiment shows that the product AB+ is formed with considerable 

vibrational excitation. Then all that one may conclude is that, for the 

reverse reaction, the use of vibrationally excited AB+ as a reactant will 

enhance the rate of formation of BC in an essentially equilibrated 

Boltzmann internal energy state distribution (characterized by the above- 

47 

mentioned temperature ^ ) . One cannot rule out the possibility that 

(for example) an increase in the translational energy (E T ,) would not be 

still more advantageous for the production of BC, if no restriction is 

placed on the degree of excitation of the BC. Classical trajectory 

calculations would be useful to explore such questions in detail for 
49 

specific systems 

The more detailed information is provided by a specific "elusive" 
experiment. Here one considers the distribution of products for a sharply 
defined initial state. The prismatic plots, as in Figs. 3 and 4, represent 
a complete summary of the energy disposal for the "forward" exoergic 
reaction, in a specific, elusive experiment fi.e. , for a narrow range 
°f fy, f R , or, more commonly, at a narrow range of E y , E R ) . A prismatic 
plot represents an enormous amount of detail, more than we would normally 
wish to know. Yet, in recording only the energy trajectory of the locus 
of the maximum in the contours in Fig. 4b, we have sacrificed too much 
detail. What is needed is to characterize succinctly the contour map 
in the region of the neighborhood of its maximum. In addition to its 
coordinates, one requires its "strength" (or magnitude) and its "sharpness" 
(or curvature). This will be considered in the next section. 
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4. CHARACTERIZATION OF A PROBABILITY DENSITY ENERGY CONTOUR MAP: 

"MOST PROBABLE STATES" 

From the discussion in the previous section it appears highly 

desirable to attempt a eencise characterization of an internal energy 

contour map, i . e. « to consolidate the overabundant detail of such a map 

J 50 

(cf. Fig. 3) by parametrizing its most significant structural features. 

Perhaps the simplest characterization is the one originally used 

by Polanyi and others 3 to summarize the chemiluminescence results. The 

data are "coarse-grained", by considering the relative rate of formation 

of product AB* in each excited vibrational state, irrespective of 

rotational excitation, and then plot the results as -k(f v t) v s^ f v i (where 

v v 13 

f yf denotes the fraction of the available energy in product vibration). 

Then one calculates the average value of the fraction in vibration, say 
?^,This implies the average fraction of the available energy released 
into relative trans lation, by difference, to be = 1 - F y i. 

In view of the trend 1,24,25 ' 48 toward the use of the more detailed 
internal energy contour maps, we should consider the concept of the set 
of three averages {7 T , F y , F R ) . For this we need to locate the "center 
of gravity" of a triangular probability contour map such as that of Fig. 3. 
This result will obviously be different than that for a right triangular 
map of relative rate constants. 1,24,48 For the latter, one could carry 
out the averaging independently in the two orthogonal directions (f y , fj^) 
to obtain T„, F D . In the absence of experiments, a classical mechanical 
Monte Carlo computer simulation (for given E T , E y , E R ) can readily yield 
number distributions in E yl and E R , (considered as independent variables) 
from which the first moments can be found, yielding Ey, and E Rf and thus 
7" and T di . (Under certain conditions, the first few moments of a p.d.f. 
can be directly obtained, with good accuracy, and with relatively few 
trajectories, without the necessity of computing the entire p.d.f. itself) . 

However, the thus -obtained average, fy, , is not the same as the 
previously mentioned average value of the fraction in vibration, calculated 
irrespective of the rotational excitation. It is also not true in general 
that Ty, + F r1 equals 1 - T v . Clearly there are many pitfalls (or at 
least ambiguities of notation) in discussing such averages, and presumably 
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there would be similar caveats applied to any "averages" over triangular 
probability contour maps,, such as Fig. 3. These considerations suggest 
that the characterization of such a map in terms of its "shape" near 
the. maximum is preferable to one in terms of average fractional energy 
disposals . 

Before proceeding further , however, it is worth noting that for 
all systems thus far studied (either by experiment or computer simulation) 
the range of relative probabilities w encompassed by the map is rather 
large, i .e. . w is usually a strongly varying function of (f^K It is 
therefore convenient to work with the logarithm rather than the original 
function itself- 

Accordingly we define a new function 1 (Ey, E^) by the relation 

w(Ey,Ej^) = ®^P [=•! (Ey ,E |P ] (5) 

where the w is the probability of Secs. 2 and 3. Similarly I may be 
defined in terms of any other pair of independent energy variables, such 
as E t , Ey, etc. 

In view of our preference to represent the results in terms of 
reduced variables, e. g . , in the triangular contour maps such as 
Fig. 3, we redefine I accordingly (noting the differential range 
relationship of Eq. 4) , so that 


1 ( £y s f jp = 

— toO)(fy S f R ) 

(6a) 

I ( f j f ifyj) = 


(6b) 


and so forth, depending upon the choice of independent variables. In what 
follows, to simplify notation, the two independent variables will be 
designated x,y. 

The locus of the "most probable final state", , y^ may then 
be determined by solving 

[9I(x,y)/9x] - . = [3I(x,y)/9y3 , * = 0. (7) 

The sharpness of the mimmura in I(x,y) is governed by the magnitude 
of the three second derivatives or "fore© constants", here designated 
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defined in the obvious manner: 


k xx 5 C3 2 ltx,y)/3X 2 ] Cm) w 

x ,y 

(© 

k xy . [ 3 2 Jfx.yy3x 3 y] m) .etc. 

x y 

In the neighborhood of the "most -probable state"; , y^ m \ 
we can expand I(x,y) and truncate beyond the quadratic terms. This 
should provide a "poor-man's characterization" of the contour map 
of the p.d.f. at the specified E. Thus we write 


I(x,y) = 


I(x 


(m) „(m) 


y ' * 7 k xxt x - x(m) 1 2 + k w tCx-x ( ' n, Vy-y ( ' m ' ) )] 


xy 1 


1. r (in),' 
* 1 


(9) 


Unfortunately, even the truncated expression (9) requires a knowledge 
of six parameters (two for the locus of the maximum, one for its strength 
and three force constants"), all of which may be expected to be E-dependent. 
One can simplify the parametrization of the I(x,y) surface by choosing a 
new pair of coordinates, say u,w, such that 1 


I(u,w) ■; 1(0,0) + J k u u 2 + j k w w 2 . ( 10) 

Such a normal mode" transformation is readily effected by diagonalizing 
the quadratic form (9); the details are briefly summarized in Appendix I. 
Since the transformation involves a rotation it requires the specification 
of one angle. Thus the three original force constants have effectively 
been replaced by two new ones plus an angle. The hope is that this angle 
will be only slightly energy- dependent, and so there may only be four 
paramters which vary significantly with E, of which two of them (the 
force constants k u , k w ) might be expected to be fairly insensitive to E 
( cf . the implications of the hypothetical Fig. 4a). 

There are, of course, intermediate levels of detail, providing 
more information than just a "poor-man’s" specification, yet not as 
detailed as the whole probability contour map. One may specify a 
distribution for one independent variable only, either by holding the 
other variable constant or by integrating over all possible values of 
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the other variable,, as described earlier. The former corresponds 
to a cut along the prismatic plot, the cut line being the line of 

constant f . In particular, such a cut can be made for f = f ^ , 
x xx 

showing the distribution with respect to the other two variables, or 
aleng m edge (£ x ® 0) , The latter is moro eommen (e£, Mg, 2), 

Note that the most probable value of f^ can be determined correctly 
from a cut at f£ m ^ but not necessarily from an edge cut (or a cut 
at any other value of f y ) . Only if the two variables are truly 
independent (!*£.. , if k = 0) or "uncoupled", are two cuts sufficient 

Ay 

to determine the distribution. 

In the next section we shall consider the implications of the 
shape of these p.'d.f. *s with respect to information content. 


' I * 1 ' » 

> 



\ 
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5. INFORMATION CONTENT OF A PRODUCT STATE DISTRIBUTION 

One notes upon inspection of a smooth distribution of product 

energies (such as the contour map of Fig. 3) or a cut through such a 

■■^.,^fc#'3itxibu.t i on, (cf t .,Fig.. 2 ), r that there is more "information content" 

in a narrow distribution (where most of the reaction probability is 

concentrated in the vicinity of the "most probable state") than in a 

broad distribution. Alternatively stated, there is less "missing 
52-55 

information" in a narrow distribution. If we know that the 

distribution. is narrow we predict with more confidence that in a future 

"experiment" most of the products will be "found" near the maximum of the 

distribution. The broader the distribution, the less our ability to make 

useful forecasts, as we are missing too much information. Clearly, a 

uniform distribution represents the maximal state of ignorance to predict 

the outcome of a given "experiment". A narrow distribution (the analog 

of a loaded die) implies that the situation is biased, favoring our ability 

to forecast the result of a future "experiment" with more confidence. 

52 

The arguments above apply, of course, to any p.d.f. ; Shannon and 
53 54 

others ’ have formalized these arguments by the itroduction of a 
quantitative measure of the uncertainty associated with any distribution. 
The more the information content of the distribution, the lesser our 
uncertainty about the outcome. In an "experiment" with a definite 
outcome, the uncertainty is zero. In an "experiment" with a- set of 
equally probable results, the uncertainty is maximal. 

The concept of the information content of a set of discrete 

52-55 

results of experiments has been fully dealt with in the literature, 
based on probabilistic and statistical considerations. It is noted that 
the number of ways in which N experiments can result in n different 

outcomes, with m. experiments yielding outcome number 1, m. experiments 

* 56 * 

yielding the i'th outcome, etc., is given by 

• _ i 

W = NI/ m.! ( 11) 

i=l 

where ? m. = N . 
i=l 1 

Let P. * m./N be the inherent probability of the i’th outcome. Then 

W = N !/” (NP.)! 

i=l 


( 12 ) 


then 
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If all outcomes are equally probable* =...=— 


W = N!/[(N/n) ! ] 


n 


(13) 


so that for large N, n we obtain (using Stirling’s approximation) 


log W = N log n. 


(14a) 


On the other hand, if only a single outcome is possible, i.e. . 

and all » 0 (k 4 j), then W = 1 and 


say P. 
} 3 


1 


log W = 0. 


(14b) 


The quantitative measure of information is usually taken to be 

I = log W (15 ) 

such that the uncertainty associated with a pair of independent "experiments" 
is the sum of, .the: independent uncertainties., , , 

Eq. -12 cars .toe rewritten for large ,N, n, and combined with Eq. 15 

to yield: 

I a — N E P. log P. . (16) 

i=l 

This is an extensive quantity. To obtain the result "per experiment", 
this is usually divided by N and written 

'•I = —IP. IgP. ( 17 ) 

X 1 

It is usual to consider the logarithm to be to the base 2, here 
designated Ig. 

Note that I is a non-negative quantity. It ranges from zero 

( cf . Eq.(14b)), when the outcome is certain, to its maximum value 
1 

(when P^ = — ) of Ig n ( cf . Eq.(i4a)) and so is a measure of the "missing " 
information content. 

In the subsequent sections we shall develop the concept of the 
information content measure which is appropriate to the different 
representations of product state distributions (both global and 
specific). Towards this eventual goal we consider in the next section 
the construct of the averaged transition, probability. 
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6. TRANSITION PROBABILITIES: DISCRETE DISTRIBUTIONS 


Let n and n ! be the sets of quantum numbers required to fully 

specify the quantal states for the reactants and products respectively. 

In general, the experimental arrangement is such that one is unable, in 

57 

principle, to fully resolve the internal states . The experimentally 

indistinguishable states are collected together in the group y . 

Suppose there are g^, final states in the group y'. The averaged 

transition probability from an initial state in the group y to a final 

25 

state in the group y' is then given by 


tt(rV) = (g y g Y .) 


-1 


E I P , 
n’ey'/iey n * 


Cl 8) 


Here P , is the detailed, quantal, state-to-state transition probability] 
the notation means a summation over those states n included in the group y . 
Subject to the reservations discussed in' Appendix II, microscopic 
reversibility ’ * holds, i.e. 


w(y^y') aai.CyW 3 = w(y,y') C 19) 

The limitations of Eq. (19) should be clearly understood. What is implied 

is that only one u)(y,y f ) is necessary to specify both the forward and the 

reverse probabilities. This does not imply that a matrix of elements o>(y,y') 

is symmetric. (In fact such a matrix will not, in general, be square). 

Nor that if u)(y,y') is small then the rate of the reverse reaction from 

state y ! to all possible states is necessarily small. 

In any particular experiment (or in any particular plot) the 

specification of the initial state may use a coarse-grained scheme 

that is less detailed than the groups of indistinguishable states y, 

(say,r s all states having internal energy in a given interval) . Let r 

be the set of possible initial groups of states. If P^ is the probability 

of finding the system in the initial group y,(yer), the transition 

59 

probability into the final group y ? is thus 

P Y » = ^ uCy')^, 


( 20 * 1 ) 
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where 


«(Y’) 


to(r,y ! ) 


(Note that P , is just the fractional or relative rate 
of reaction into the product group y f ). When the initial state is 



definitely in a particular group y 

-1 


(I*®.- P v s 1) 


P v t = w(Y»Y' 

Y 


= g„ 


E E 

n'ey^ney 


Y 

P , 
n’ ,n 


( 21 ) 


corresponding to the usual procedure of averaging over initial states 

and summing over final states. 

One can consider either of two (tirae-reversal invariant) 

probabilities,, either w (y^Y 1 ) * the averaged probability (which, as is 

59 

clear from Eq. (18 ) is bounded by 1), or the yield 


Y(Y.Y') ■ g^i^YaY') « (22) 

Y Y 

When the initial state is prepared without any attempt to resolve the 

states within the group 'T then P is sicply proportional to g and 

Y Y 

hence 


p v , = s'(p v /*J Y'(; r , r’) (23) 

T y T Y 

(Note that ec is the equilibrium condition for the microcanonical 

ensemble, i.e. for an ensemble where the total energy is in the range 

61 62 

E to E + dE, all quantum states are equally likely. * ) 

Hie formal theory supporting these results is summarized in 
Appendix II. 

As an example of such averaging procedures, consider global data. 
We restrict attention in this section to the discrete case. ' A particular 
entry in . a poor-man's P-natrix is 




Using Eqs. (21»2S} 9 this expression becosses; , 


(24) 
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Here the primed summations are over y and y' (or n and n') in the range 
corresponding to the entry in question. (The analogous results for the 
case of continuous internal energies are developed explicitly in Sec. 9). 

Next we consider a specific-type experiment, where the initial 
conditions are identical to those specified previously in u> but the final 
distribution of products is known. The triangle plot (cf. Fig. 3) can then 
be expressed in terms of 

u(y') = E' P iii(Y)Y') ,yer . -> (26) 

Y r 

Thus 

«o(r,r*) = E’g^Y')/ E/g Y , , y'er* . (27) 

y' t' 

All the information down a given vertical slab of a bar plot 
(whidi corresponds to specified initial conditions) can be obtained from 
a triangle plot for the appropriate specific- type experiment. The necessary 
modifications for continuous distributions are given in Sec. 9. 
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7. INFORMATION CONTENT : SPECIFIC EXPERIMENTS 


In this section we consider the information content of the product 
state distribution for the case when a partial resolution of quantum 
states is possible. (Continuous distributions are deferred until See. 9) , 

The averaged (over the group y') probability of finding the products 
in the state n* is (cf. Eq.(20)) 


P , = P ,/g , 
n ’ y ’ 6 y ’ 


(28) 


For the case of an "inclusive" specific experiment this becomes 
(cf. (20)) 

P n , = E’P w(y,Y*) = to(y') 

Y T 

and for a non-inclusive ("elusive) specific experiment 


(29) 


P n , ■ w(y,y') 


( 30 ) 


In both cases 




(31) 


It should be clearly stated that (28) assigns equal probability to all 

the states within the group y'. (Recall that final states within the 

group y' are, in principle, indistinguishable, under the experimental 
57 

conditions used .) . 

The information associated with this product distribution is, 
from Eq. (17) 


E P .IgP 


n' 


n’ 


= -E 2 (P v ,/g v ,)lg(P v i/g Y .) 
y' n’ey' Y Y Y Y 


(32) 








21 


where in the second line the inner summation is over g , identical 

Y 

terms. The fourth line serves to define a quantity Ip(y'), to be 
referred to shortly. 

Given the experimental product state distribution and thus 
the P /, one can obtain I p via Eq. (32) as 

I p = - <lg(P Y ,-/g y .)> (33) 

where the average is over the final product distribution. Ip is 
clearly bounded: 


0 s? I p s? 


lg(£ g v i) 
Y 


(34) 


The lower bound corresponds to all the products being in one particular 
group while the upper bound is the microcanonical equilibrium distribution, 
i.e. when P ,« g^ t (so that the probability of the products being in 
the group is proportional to the number of states in the group 
( cf . Sec. 6)). These results are essentially those already given in 
Sec. 5, Eqs. (14a,b) with Eq. (17) for I. 

i i ■ i 

Ip is a measure of the information content of the whole distribution 

We can also define an information measure of a particular outcome, say y'. 

62 

This quantity has been termed the "surprisal" . It has been anticipated 
in Eq. (32) and is defined by 


ip(Y') - -ig(P Y ,/g Y ,) » -IgaCy') . (35) 

Thus Ip is the average value of the surprisal: 

Ip = E P v> I P (Y ’ ) = V 1 gt“('r')] = < I p(Y , )> • (36) 

Y» » y ^ * 

Eq. (36) and its implications noted below are the main results 
of this section. The information content of a distribution is the 
average of the information contents of the particular groups, each 
group corresponding to a, possible outcome of the reaction. (By 
definition, we are unable to distinguish between the members of any 
specific group) . 
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It is of interest to compare the surprisal for the two types of 
specific experiments, the "elusive'-’ and the "inclusive". For the former 
the surprisal for a particular outcome is 


I p Cy sY * ) = -“lgD»>(Y ,Y ')] . (37) 

Here I p (yjY*) is the elusive information measure of the outcome y ' for the 
initial group of states y. In an inclusive experiment one necessarily 
first averages over the entire distribution of unresolved initial 
groups of states 


w( Y ') - S # P w(y,y’) (38) 

Y T 

and then obtains the information measure ( cf . Eq. (20)) 

Ip(Y s ) = -lg[w(Y')3 = -lg[S ff P »(Y*Y')] (39) 

Y T 

It follows from the inequality for convex functions 54 ’^ 3 that 


I P (Y 5 ) P Y Ip (y,y’) . 

•v/ ' 


(40) 


If we now average these surprisals over all final states to 
obtain the information associated with the overall product distribution, 
following the procedure of Eqs. (32) and (36), we see that the value of 
I p for the inclusive case (obtained by averaging the £.h.s. of Eq. (40) 
over all yj) is greater than the I p for the elusive case (from the r.h.s, 
of Eq. (40)), The latter is closer to being the "true" or intrinsic 

: i 

information associated with the product distribution for the specific 
experiment. 

Instead of considering the product information measure 1 , one 

can consider the information defect with respect to a microcanonical 


equilibrium situation f e. g . . Eq. (33) vs_, the upper bound of Eq. (34).). 

We propose to measure the information defect by the average 
deviation of Ip(Y') from its microcanonical equilibrium value 


&I p = <[I p (y')] - X p (Y')> (41) 

We note that for a microcanonical equilibrium situation at a given 
total energy E in the small interval between E and E + dE, all quantum 
states are equally probable (a microcanonical ensemble), hence 

(P Y 1 3 gq “ Syl/^Sy* 


(42) 
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and 

t I pCY') ] e q - "IgtS g y ,) = -lg[tP Y ,)gq/g Y ,l t«) 

Thus, the information defect is 


Up ■ P Y , lg CP Y ./CP Y ,D eq J . (44) 

It is expected that Al p >, 0, with equality holding only for equilibrium. 

This follows from Shannon's Lemma: If E P , = 1 and E,(P ,) ori = 1, then, 

_ 1 y • Y y Y ec l 

since lg x ^ 1 — x" , with equality if x = 1 , then 


Al_ *E P ,[1-(P ,) /P , ] = E P . — E (P , ) =0 (45) 

P y' Y Y eq y' y . Y' y , Y 1 eq 


Thus, for any non-equilibrium distribution there is a net positive 
information defect, relative to the (microcanonical) equilibrium 
distribution. ( A small AI^ Implies a more uniform distribution, a large 

Alp a more sharply peaked distribution.) 

Based on the results following Eq. (40), it is seen that the 
information defect for an "imperfect", inclusive type experiment 
will be less than the intrinsic value of the information defect as 
obtained from an ideal elusive type specific experiment i.e. . 


^ A *P^ elusive ^ 


(AT). > 0: 

v P' inclusive 


(46) 


Of course, (AX ) .. . is the more correct characterization of the 

p elusive 

specific experiment.” It is bounded from above by the information defect 
for the (microcanonical) equilibrium distribution. 

The equivalence between I and the thermodynamic concept of entropy 

53 P 54 59 

was shown by Jaynes and Khinchin, based essentially on Shannon’s 

work. This will be employed and extended to deal with the present 

application in Sec. 10. 


8. GLOBAL AND SPECIFIC REPRESENTATIONS, 


We have introduced a hierarchy of representations, depending on 
the resolution obtained, which can be summarized as follows. 

(a) Global elusive . This is the most detailed characterization 
of the reaction possible (in principle) at a given total energy and 
specified experimental conditions. It corresponds to a complete 
specification of u>(y,y 8 ) for all possible groups of inital and final 
states. w(y ) is the "rich experimental^ tV°P«=> matrix. In contrast, 
the rich quantum theorist considers P , as the elements of the P-matrix. 

ill ^ Jfi 

However, under given conditions the experimentalist can only determine 
the average of over ney and n'ey 0 (cf. (18)). 

(b) Specific elusive . The most detailed summary of the specific 
ideal experiment , showing the distribution of products oj(y,y ! ) for some 
initial group of states y. A specific elusive experiment is equivalent 
to a column of the global elusive representation. 

When one goes over to a continuous distribution of products 
the classical limit) the specific elusive representation is the 
result of a classical trajectory calculation. Figures 3 and 4 are 
specific elusive representations. 

The current state of the experimental art has not yet reached 
the stage where both the initial and final states can be fully specified. 
The following three representations take cognizance of this fact, by 
averaging over either the detailed distribution of the reactants or of 
the products or both. We can thus consider the following. 

(c) Specific inclusive . The detailed distribution a>(r,Y') 
of products from some initial ensemble V° Thist is the distribution most 
commonly obtained from chemiluminescence experiments and, as discussed 
before ( cf . (46)), provides an information deficiency that is nearer 
to the equilibrium value than that obtained from elusive experiments. 

Very recent results 4 have now provided an insight into the dependence 
of w(r,y s ) on the initial ensemble P. Alternatively, one can also consider 
“(Y,r*) , the influence of the reactants' state selection on the probability 
of formation of products . ■ 
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(d) Global inclusive . A poor man’s P-matrix, in the form of to(T,r'), 
say w(f,p f^,), as shown in Figs. 1 and 5. An entry in the global inclusive 
representation can be obtained by an averaging of o)(r,Y') (or ujCy > T ’ ) ) over 
Y',Y’er’ (or over yeT) > as shown in Fig. 5. Along the same lines as the 
proof of (40) one can prove that the surprisal associated with co(r,r') 

f i.e. , — lg o)(r,r') ) is larger than the average of the surprisals of u>(r,Y')- 

Thus, an equivalent statement to (46) obtains. A poor man's E-matrix 

is nearer to an "equilibrium" distribution than a rich man's P-matrix. 

Any averaging necessarily reduces the information deficiency and brings 

the resulting distribution nearer to the equilibrium (or statistical) limit. 

2 8 29 

To a certain extent, the success of the statistical theory ’ depends on 
this degradation of information by averaging. 

(e) Inclusive , . The direct experimental determination of a column 

of the global inclusive P-matrix. Here one is concerned with transitions 

between some averaged set of reactants (T) to an averaged set of products 

(T'), reporting the dependence of u)(r,T') either on r' (the averaged 

distribution of products) or on r (the averaged role of reactant excitation). 

7-9 

Molecular beam experiments with velocity analysis provide the fy, 
dependence of w(r,f T ,). The pioneering studies of the influence of reactant 
internal energy were also of this type. 


IIV).',*; 


i i-j ) . 
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9. INFORMATION CONTENT OF A CONTINUOUS DISTRIBUTION 


As discussed in the earlier sections, individual quantum states are 

often not resolved and one observes continuous rather than discrete product 

state distributions. Problems involving continuous distributions are 

53 55 

reputed to be more difficult * because the concept of "equal probability" 
is less clear cut. In the present application, we have dealt with this 
problem already even for the discrete case. (Recall that the assumption 
expressed by Eq. (28) is that we assign equal intrinsic probabilities to 
each state in a group which is, in principle, unresolved experimentally.) . 
Extension to "valid" continuous distributions requires nothing more, 
beside the conversion of summations over quantum numbers to integrations 
over energies. Of course, proper care is required to identify equi-probsbl i 
quantum states; thus appropriate density of states factors are needed. 

Our goal is to obtain an explicit expression for the average value 
of the suprisal. Ip, for the case of a continuous product (or reactant) 
state distribution. We shall consider E£ as the (continuous) energy corresponding 
to the formerly discrete variable, say £, which is now (in the "classical case") 
continuous. Typically this would be the translational energy E^* ( cf . Sec.' 2). 
The group f is the group of states with energies in the interval E^ to 
E^ + dE^ and g^ = ^(E^ is the density of states, so that g (E^) dE^ is 
the number of states in that energy interval. = P(E^) is also now a 
continuous p.d.f. , with P(Ep dE^ being the (averaged) probability that 
the product energy lies in the interval to E^ + dE^. (Similar considerations 
apply to reactant state distributions.) ,, 

The surprisal I^ff) = IpCE^) is also now a continuous function of E^, 
with the above interval having the information measure IpfE^) dE^ . By 
analogy with Eq. (35) , we have 


Ip(E f ) * -lg[P(E f )/g(E f )] 
= — lg[io(E f )] . 


(47) 


Here ^(E^) is the averaged transition probability to the group £, where 
the discrete averaging of Sec. 6 ( cf . Eqs. (18) and (20)) is now an average' 
over a small energy interval. The transition probability to the group f is 



given by 
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P(E f )dE f = co(E f )g(E f )dE f C48) 

(Of course, in practice, this is how one obtains to(E^) from the observations.) 
Appropriate densityr-of-states factors (for di atomics in the rigid. rotor- . 
harmonic oscillator approximation) suitable for either bar or prism plots, 
have been provided by Kinsey. ^ 

The information content of the distribution is, as in Sec. 7, the 
average of the surprisal (cf. (32)) 

I p = <I p (E f )> = /dfu)(E f )I p (E f ) 

= /dE f a)(E f ) g(E f ) I p (E f ) (49) 

= — /dE f P(E f ) lg[u>(E f ) ] . 

In the case of more than one energy variable, e. g. . Ey, , Ep, as well 
as E,p t , extension of the above should be readily accomplished, taking 
advantage of the approach of Sec. 4. 

In practice the continuous integration in Eq. (49) would be carried out 
as a summation over energy ,r bins": 

I = -E P(E )[lgu,(E ,)] .AE . (50) 

AEy' 1 

(The results should become , asymptotically at least, independent of bin size). 
Global results imply a knowledge of the dependence of the reaction 
probability upon both initial and final energy variables. In the simple 
case of a bar plot (cf. Sec. 2) one can use either the averaged transition 
probability o)(E T ,E^,) or the yield Y(E^.,E T ,). From Eq. (23) . 

P(E t ,) = /dE T f(E T )y(E T ,E T ,) , (51) 

where g(E T )f(E T )dE T is the fraction of initial states with E T in the range 
E t to E t + dE T . Similarly from Eq. (20) 

wtE T ,) = /dE T g(E T )f(E T )a)(E T ,E T# ) (52) 

Corresponding results hold for the integration over E T , . 


• ; ' V . . J 



10. ENTROPY DEFICIENCY OF A PRODUCT STATE DISTRIBUTION 

There is a well-established equivalence between information content 
and the thermodynamic concept of entropy . 52 " 55 j n the case of 
equilibrium distributions the situation is clear cut, but in the present 
non-equilibrium context it must be admitted that our recipe for assigning 
probabilities by the equilibrium rule "all quantum states of the same energy 
are equally probable" is less secure. Nevertheless this is the simplest 
means of extending the more familiar equivalence to the problem of the 
characterization of a non-equilibrium product state distribution . 65 

Assuming the validity of the equivalence, we next establish the 
practical connection between Ip and S in terms of numerical magnitudes 
and units. From Eq. (17), if I, the "missing information content", is 
based on the logarithm to the base 2 , jLe^ I a ig w, then to obtain the 
usual result for the entropy, Slk log t'J from Eq. (15) one must write 

S = (k log e 2) I p (53) 

A measure of the specifieify of tit. reaction is the entropy 
deficiency in the product state distribution, which can be expressed 

AS"= (k logg2) Alp (54) 

with given by Eq„ (44). In practical terms, (i.e . a for one mole of 

reaction), the entropy deficiency is thus AS'(e.u.) ° X.38AI . We recall 

66 P 

that the units of information are bits, with 1 bit = Zg2 . Since the size 

of the entropy deficiency AS' is perforce limited (i.e., the bound being 

commensurate with the magnitude of standard entropy change of reaction AS°) 

there appears to be an intrinsic upper bound on AI (bits). • ■ 

P 

The entropy deficiency would be zero if the product state 
distribution were that predicted by a phase-space or statistical (density- 
of-states) theory. For all actual situations AS ' 1 is intrinsically positive. 

In principle, the product state distribution of a elusive specific experiment 
would be suitable for analysis to ascertain the entropy deficiency AS' 
of the reaction (although, due to imperfect reactant-state selection and 
product resolution, it would yield only a lower- limit estimate of AS). 

0 ns can recognize a hierarchy of experiments leading to a range of 
AS* values , In the absence of data, forced to use a statistical theory, we 
start with AS ° 0 . As we continuously sharpen the experiments the computed 
AS will be larger 5 finally the limit ©£ well-defined reactant and .product 

states (from a global elusive 
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representation) would yield the proper intrinsic value of the entropy 
deficiency for the reaction at the total energy E. Presumably one would 
then be in a position to consider the overall dependence on E. 

Clearly the concept of an entropy deficiency AS* should also 

68 

be a useful one in characterizing the poor-man's P-matrix. Whether it can 
be illuminating physically is considered in the second paper in this 
series. Its relationship to the heuristically valuable concept of an 
entropy of activation for a canonical system at a temperature T remains 
to be established. 


i i * 
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11 . SUMMARIZING REMARKS 

This paper has attempted to characterize the product (and reactant) 
energy distributions obtained in reactive scattering experiments or computer 
simulations . Suitable forms of plotting results of global or specific-type 
experiments have been suggested in Secs. 2 and 3 (see Figs. 1 , 3 , 4 ). Means 
of parametrizaticn of specific type results in terms of an expansion around 
the most-probable states were developed in Sec. 4, with the most compact 
form being that of Eq. (10). 

The question of the information content of a discrete product state 
distribution was discussed in Secs. 5-7 with the important results being 
those of Eq. (17) for I, Eq. (32) for Ip and Eq. (44) for Alp (the 
information defect). Sec, 8 summarizes the types of representations used 
and their interrelations and stresses the concept of the degradation of 
information by averaging. The information defect, AI p , is shown to decrease 
with any additional averaging. The intrinsic Alp can only be obtained fro® 
elusive experiments. Any less detailed results yield only a lower bound. 

In Sec'. 9 the generalization was made to continuous distributions, with 
the results for I p given in Eqs. (49) - (50). 

Finally, the concept of entropy deficiency AS^was proposed in 
Sec. 10 and quantified in Eq. (54). 

Obviously there are many theoretical as well as practical 
ramifications of the present, rather preliminary, investigation. These 
are currently being pursued. 
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APPENDIX I 


There are at least three ways in which one can plot a function 
I (E^,, Ey, E R ) whose variables satisfy a conservation relation, i. e. ? 

Eq. (2). (a) An equilateral triangular plot, as discussed in Sec; 3. 

(b) A cartesian axis plot, as used in Refs. 24, 48, etc. and (c) a 
triangular plot where the skewing angle is defined so as to diagonalize 
the function I about the minimum. 

The advantage of choices (b) and (c) is the ease of diagonalizing 
the function. For the cartesian case, one requires a unitary matrix U, say 


f zosq sine \ 

U = . 

\-sme cose / 


’ e sine 
•sine cose 

Here 6 is determined such that 


(A.l) 


U 



(A. 2) 


In this case, the maximum is determined by its two coordinates, by the 
angle 6 and the two diagonal derivatives k and k . Case (a) is analogous 

W 

to case (b) but requires a pretransformation from the skewed coordinates 
to cartesian coordinates. 

Case (c) is analogous to the method of plotting the potential 
energy surface for collinear collisions, with the exception that here one 
chooses the axis not to diagonalize the kinetic energy but what corresponds 
to the potential energy. In fact, we have loosely referred to the second 
derivatives of I as the "force constants". This is not simply because of 
the mathematical analogy but because of "thermodynamic" reasons (to be 
discussed elsewhere) . 

Choosing Ey and E R as the independent variables we scale the E R 
axis by c f i. e. . to be plotted is E R /c) and tilt it by the angle 0 from 
the (vertical) Y axis . Thus 


E R = c Y, sees 


Ey = X _ Y tan0 


(A. 3) 
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where Ey is parallel to the X axis. The angle 6 is selected to eliminate 
any cross derivative, i.e. , 

= kxjr / kyy CA.4) 


and the scale c is defined so that ky = k% or 


c 2 = k /k 

xx yy 

(A. 5) 

sin 2 e = k^/k^ 

(A. 6) 

With these choices. 

I(X,Y) = I(XW Y^) + | p 2 k 2 

(A. 7) 

where 

2 >• (m).2 a r (m).2 

P = CX " X ) + C.Y - Y ) 

(A. 8) 


2 

and k ■ k^. Superscript ( m ) denotes the most probable value. 

In this way the maximum has five parameters (Ey 11 ^ , , c, 0 and k) 

but only one second derivative. 
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APPENDIX II 


We are concerned with transitions between two sets of states Psfrl 
and P — ^ . Let p(P) and p(P*) be the appropriate density matrices 

for the two sets. One can then define the rate constants for the forward 
and reverse reaction* 39 ' 60 


fe(r-> r')=Tk[p (P) S p (P 0 S+] (\xT v (p(r) } 


(B. 1) 


and 


(r r)-Trfp(r')sp(p) sf } / kTjpCr 0\ <*• 


2) 


In both cases we average over the initial ensemble and sum over the final 

ensemble. When the states in the group S/ cannot be resolved under the 

° 59 

given experimental conditions they are taken to be equiprobable. The 
yield function and averaged transition probability are defined by 


y(r->r')-Tv-{p(r)sprp')s^ ] 


(B.3) 


and 


^ (r->rO = V 0v> 

The yield function^ P— >P!) represents the sum of the transition 
probabilities while the probability function t)(P— >r # ) is the averaged 
transition probability, evaluated at a specified total energy E. 

If p( P ) and p(P*) are separately invariant under time reversal, 
( not necessarily always the case) then both Y(P*^yP ) andfi^CP-^p^ 
are also invariant, i.e. . 

y(r->r') = y(rUr)=Y(r\r‘) (B . 5) 

and similarly for t^Cr^P^). When p(P) is a microcanonical distribution 
one can use the shorthand notation Y(E) (or /^(E)). 

In general, it has been shown ( i.e . . Ref. 59, p* 143) that the 
time-reversed form of an operator, here denoted by a bar, is given by 


p(p) = 9pf(r')© -< 


(B.6) 


where Q is the time-reversal operator. Since the density operator 
is Hermitian, it is time reversal invariant if 
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p(r)=©pCr)© '=p(r), cb.t, 

Using the properties of the trace one can go through the following sequence 
of operations and show that Y(P-~>P^) 21 ( P O • 

if p = F . 

T r {p lr ) S p O'" ') S+ } - Ty {p(n' ') <, + p (P ) S } 

(B. 8) 

-"k{ep(rO 0 H s^0p(n) (r)e Vep(r)e H :sp} 


From Eq. (B„6) (or Ref. 59, p. 147), 

s = es + e r1 


and since 



<S=9’k+©. 


(8.9) 


(B . 10) 


Thus the last term <^n the right hand side of Eq. (B.8) becomes 

T^fp(r ( )Sp(P)s+y 

seen to be the same as the first term (the l.h.s.) with P* and P ^ 
interchanged. 

y (p _>]->') -Ty tp(r)sp’{r') 

- T v {p(r 0 s p ( r) ^ p ) CB ' n) 

as expected, so 

y(r-> r 0 = V T r)= y C v > r ')=Y^P r ). (b;5) 
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Legend for Figures 

1. (a) Contour map of w (jj^E t ) at a fixed value of the total energy, 

E = 10 units, for a (hypothetical) exbergic reaction for which **AE q *» 5 

units. Thus E^ and E^,, range “from 0 to 5 and 0 to 10 units, respectively. 

(b) A square-faced bar plot w(f T »f T ») for the same reaction, showing the 
dependence of the contour map upon the total energy, over the range 

10 E < 50 units. The "trajectory 1 * of the maximum of w as a function 
of E appears as the heavy curve passing down the length of the bar. 

For the purpose of drawing the results in (a) and (b) they were scaled such 
that, at the maximum, a) has been assigned the value 1, i.e . , (o' 1 -(E) = 1* 

(A similar comment applies to Figs. 3-5.) An actual experiment (or compu- 
tation) would also provide the energy dependence of the peak value, (E) . 
Thus, besides Fig. 1(c) one would also require a figure showing vs. E . 

• / V / V 

(c) The coordinates of the loci of the maximum, i.e ., and f^, , plotted 

vs. total energy E over the same range as in (b) . 

2. Comparison of the effect of E^ with that of E T upon the total reaction 

+ + 

cross section for the reaction ^(v) +0 ->- NO + N (AE q = l.leV). The 

Ey dependence is taken from Ref. 17, while the E^, dependence is from Ref. 43. 

3. Triangular product contour map of w(f T ,, f v ,, f R , ; f v , f R ) for some fixed 
(but here unspecified) f^ , f R , E. In the example shown, the locus of the 
most probable w is f^,, = 0.15, f R , = 0.25, fy» = 0.60 (most of the avail- 
able energy going into product vibrational excitation). Here and in Fig. 4 
the triangular maps represent the ’results of specific elusive experiments. 

(cf . Sec. 8). 

4. (a) Prismatic plot of w(f^,,, f^,, f R , , fy, f R ) showing the dependence of 

the triangular contour map of Fig. 2 upon the total energy E (5 £ E £ 20 units), 
(b) The coordinates of the loci of the maximum, i.e. f^°) fffl, * 

1 V R 

plotted vs. E over the same range as in (a). 
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5. ’’Decomposition” of a single point on the square contour map wCf^, f^,) 
into its two ’’components”, i.e. . cuts through the two triangular con- 
tour maps U)(f v , f R , a; f T ») and wCfy** f R * s b; f T > . The global maps 
here and in Fig. 1 are global Indus ive representations. Thus these 
triangular maps represent specific inclusive experiments. 
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